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ABSTRACT 

Context. Based on current theoretical models, giant planets can form via gravitational instability or core accretion. The 
terrestrial planets are "failed cores" in the core accretion paradigm that do not evolve into gas giants. Planets residing 
in binary systems put a strong constraint on any planet formation theory as the model must work in binary systems 
too, in spite of the strong perturbations from the secondary star. 

Aims. We examine the first phase of the core accretion model, namely the dust growth/fragmentation in binary systems. 
In our model, a gas and dust disk is present around the primary star and is perturbed by the secondary. We study 
the effects of a secondary with/without eccentricity on the dust population to determine what sizes the aggregates can 
reach and how that compares to the dust population in disks around single stars. 

Methods. The equation of motion of dust aggregates including gas drag and gravitational forces from both stars, as 
well as the hydrodynamical equations for the gas disk are solved. We determine the velocity of these particles relative 
to the gas (or relative to micron sized, well coupled dust particles) and we construct a collision model with growth and 
erosion. Particles below the critical fragmentation velocity increase in mass by sweeping up the small dust, particles 
above this critical speed are eroded due to collisions with the small dust population. We use this model to determine the 
sizes that aggregates can reach in the disk and we also determine the equilibrium mass and stopping time distributions 
between growth and erosion. 

Results. We find that the secondary star has two effects on the dust population. 1.) The disk is truncated due to the 
presence of the secondary star and the maximum mass of the particles is decreased in the lowered gas densities. This 
effect is dominant in the outer disk. 2.) The perturbation of the secondary pumps up the eccentricity of the gas disk, 
which in turn increases the relative velocity between the dust and the gas. Therefore the maximum particle sizes are 
further decreased. The second effect of the secondary influences the entire disk. Coagulation is efficiently reduced even 
at the very inner parts of the disk. The average mass of the particles is reduced by four orders of magnitude (as a 
consequence, the stopping time is reduced by one order of magnitude) in disks around binary systems compared to dust 
in disks around single stars. 

Key words. (Stars:) binaries: general, Planets and satellites: formation, Protoplanetary disks, Methods: numerical 



1. Introduction 

Among all the planets known around main sequence stars, 
roughly 10% resides in binary systems. Al l of the s e sys - 
tems have a so called S-type configuration (jDvorakL I1986D 
where the planet orbits around one of the stars and 
the companion acts as a perturber. Most of these bina- 
ries have separations larger than 100-300 AU in which 
case the secondaries play a more limited role in the for- 
mati on, dynamical evolution and migra t ion o f the plan- 
ets (|Desidera k Barbieril 120071: iDuchenel l2010l ). However, 
there are some systems with small separatio ns (~ 20 AU) . 
These are 7 Cep with a s epara tion of 18.5 AU (Ha tzes et all 



120031: iNeuhauser et all I2007D. G L 83 with a separation of 



18.4 AU (lLagrange et all 120061), H P 41004 with a separa- 
tion of 23 AU (jZucker et all l2003h and HP 196885 with a 
separation of 17 AU (jCorreia et all [2008). The mere exis- 
tence of these systems puts a very strong constraint on any 
planet formation theory as the models have to be able to 
produce planets in these close binaries as well. 

The core accretion paradigm of planet formation 
(|Mizund . Il980l: iPollack et allll996D has three stages, start- 



ing with the coagulation of sub-micron si zed grains which 
leads to the formation of planetesimals dSafronov 119691 : 
iWeidenschilling k Cuzzl fl99l iBlum k Wurml . 120081) . The 
next stage is the fo rmation of the protoplanetary cores from 
the planetesimals (1 Weidenschillina. Il980t iNakaeawa et all 
ll983tlWetherulll990tlTanaka et all 120051 ) . Finally, the last 
stage of planet formation is gas accretion onto these cores, 
or - in case the gas is absent - gravitational encounters be- 
tween these protoplanets, which result in a chaotic impact 
phase, until o r bital stability is achieved (ICnambersl . [2001; 
iKokubo et all . 120061: iThommes etafl I2008D . 

This picture is already very complex and has many 
unresolved problems, we only mention here a few. The 
dust coagulation process, the first stage of planet for- 
matio n has to overcome the charging barri er (jOkuzuml 
2009), bouncing barrier rtZsom et al. . l2010f ). radial drift 
barr ier jW cidcnsc hillingl. I1977D and the fragmentation bar- 
rier (jBlum k Wurml . 120081 ). Pue to these uncertainties, the 
i nitial planetesimal size distribution is poorly constrained 
(|Morbidelli et all I2009D . therefore the second and third 
stages are also affected by these barriers. The presence of a 
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close companion further complicates this processes in many 
ways. 

The tidal torques of the companion generate strong 
spiral arms in the disk around the primary. Angular mo- 
mentum is transferred to the binary orbit which will trun- 
cate and restructure the disk (lArtymowicz fc Lubowi . [l99l 



lArmitaee et~all 119991: iKlev et al.L l 2008h . This dynamical ef- 
fect of the secondary has several consequences which might 
influence planet formation: it decreases the lifetime of the 
disk, increases the temperature of the disk and modifies the 
stabile orbits around the primary. 

The evolution of planetesimals is also influenced in a bi- 
nary system. The perturbation of the secondary increases 
the relative velocity of the planetesimals and/or creates un- 
stable regions where the planetary building blocks cannot 
maint a in a stable orbit as shown by e.g., |H cppenh eimerl 
(11978ft . IWhitmire et al.1 (|1998ll . lThebault et all (12004ft . The 
increased relative velocities can t hen lead to the disrup - 
tion of the planetesimals. However. lMarzari fc Scholll (l2000h 
showed that the combined effects of the gravitational per- 
turbation and the gas drag may increase the efficiency of 
their accretion by reducing their relative velocity and pro- 
duce, later on, terrestrial planets. iMarzari et al.l (|2009h cal- 
culated the relative velocity of planetesimals in highly in- 
clined systems and concluded that planet formation ap- 
pears possible for inclinations as high as 10°, if the sep- 
aration between the stars is larger than 70 AU. The re- 
gion where planetesimals can accumulate into protoplanets 
shrinks consistently for lower binary separations. 

Numerical simulations have shown that it might be 
possible to form giant planets via gravitation al ins t abilit y 
which is the other p lanet formation paradigm (jBossl ((2006) , 
iMaver et al.l (|2007l ) and references therein) , although some 
of these studies had to assume larger than 20 AU separa- 
tions (50-100 AU), or rapid cooling of the gas. We, however, 
concentrate on the initial stages of the core accretion model, 
thus we study the formation of terrestrial planets and/or 
giant planets formed by core accretion. 

Several authors studied planet formation and 



migration in the 7 Cephei system (Dvorak et al., 


2004: Thcbault et all 12004 HaghiehipourL 


2006: 


Verrier & Evansl. 12006 


IKlev & Nelson. 


2008: 


Jane-Condell et al. 12008 


: iPaardckooper et al.l. 


2008: 


Paardekoooer & Leinhardt 


. |2010D. We also consider the 



as a starting point in our study. We focus on the motion 
and relative velocity of dust particles, thus investigate the 
first stage of planet formation in this environment. First, 
we produce a surface density profile of the gas, which is in 
quasi-equilibrium and truncated by the secondary. Once 
the equilibrium configuration is reached, we follow the 
motion of 10 5 large dust particles. The small, sub-micron 
sized particles are well coupled to the gas, therefore their 
relative velocity compared to the gas is zero. The relative 
velocity of larger particles (pebbles, boulders) can be 
therefore considered as the collision speed between a well 
coupled sub-micron sized particle and the large particle 
in question. Based on this recognition, we construct a 
simple collision model in which the larger particles lose 
mass, if their relative velocity is larger than the critical 
fragmentation velocity, and increase mass below this 
critical speed. The critical fragmentation velocity can be 
determined by laboratory e xperiments and its v alue is 
usually 1 m/s for silicates (jBlum fc Wurml . I2008h . This 



way we determine the maximum particle size reachable in 
these disk configurations. We examine the properties of 
these particles to find out how far can coagulation proceed 
in such an environment, what is the size of the building 
blocks of planetesimals in a binary environment and how 
does this compare to the sizes reachable in disks around 
single stars. 

The paper is organized as follows: in Sect. 2 we describe 
the model set-up of the gas disks and the numerical method 
to integrate the particle motion and mass; in Sect. 3 we 
present our results. Finally, we discuss the results and draw 
conclusions in Sects. 4 and 5. 

2. Numerical modeling 

2.1. Gas 

We assume that the primary star is surrounded by a flat, 
non-flaring and non-selfgravitating disk which is coplanar 
with the orbit of the secondary. Therefore, we perform 
2D hydrodyna mical simulations using the FARGO code 
(jMassetl . l2000f ) , which solves the vertically integrated hy- 
drodynamical equations. This approximation is valid if the 
pressure scale height of the disk (H p ) is small compared to 
the radial distance from the star (h = H p /r <C 1, where 
h is the aspect ratio of the disk). Th e intrinsic turbulence 
of th e gas is assumed to be of a-type (jShakura fc SunvaevL 
Il973l ). the turbulent kinematic viscosity (v) is 

v = ac s H p , (1) 

where c s is the isothermal sound speed. The value of the 
a parameter reflects the strengths of the turbulence in the 
disk, typical values are between 10 -5 — 10~ 2 . Finally, we 
assume that the disk is locally isothermal, the sound speed 
is 

c s = v Kcp h, (2) 

where WKcp is the Kepler velocity. FARGO uses spherical 
coordinates (r, ip) and the origin of the coordinate system 
lies in the center of the primary star at r = 0, tp — 0. The 
coordinate system is non-rotating. 

We perform 8 simulations using different initial con- 
ditions. Our 'fiducial' simulation uses the stellar parame- 
ters of the 7 Cep system, namely we choose the mass ra- 
tio between the secondary and the primary to be 0.286 
(Mi = 1.4M , M 2 = 0.4M Q ), the eccentricity of the 
secondary to be 0.4, and their separation abin = 20 AU, 
which translates i nto a n orbital period of P = 66.7 yrs 
(jNeuhauser et all l2007f) . We list the parameters of all the 
simulations in Table [TJ We perform a simulation without a 
secondary ('single') assuming a simple power-law disk pro- 
file and determine the maximum particle size in this disk. 
Comparing the outcome of the 'single' simulation to the 
'no_ecc' simulation, in which the companion has zero eccen- 
tricity, enables us to determine the effects of the presence of 
a secondary on the dust population. Similarly, comparing 
the 'no_ecc' results to the outcome of the 'fiducial' simula- 
tion makes it possible to determine the effects of an eccen- 
tric companion on the dust population, etc. For the rest of 
the simulations, the initial parameters are chosen such that 
only one parameter has a different value than our 'fiducial' 
simulation except for the 'alpha' simulation. We found that 
with a non-reflecting boundary condition with a — 10 -3 , 
the eccentricity of the gas is not converging within 200 sec- 
ondary orbits. We found that the average eccentricity of 
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the gas disk is 0.003, 0.07, 0.17, 0.23 and 0.28 at t=0, 20, 
60, 140, 200 orbits respectively. As these simulations are 
computationally expensive (on 8 cores it takes roughly a 
week to simulate 200 secondary orbits), we decided to use 
an open boundary condition where the eccentricity of the 
gas converges well within 200 orbits. 

Mapping the parameter space in this way enables us to 
determine which parameter has a strong influence on the 
maximum particle size. 

2.1.1. Initial conditions 

We use 240 logarithmically spaced radial bins and 504 
azimuthal bins to simulate the disk, resulting in nearly 
quadratic gridcells. It extends from r; m = 0.4 AU until r out 
— 8 AU. In the case of the 'r_in' simulation, the radial bin 
goes from 0.1 AU until 8 AU and we use 351 logarithmically 
spaced radial bins to keep the size-ratio of the grids con- 
stant. The initial surface density profile is S = So r ~°' 5 j the 
surface density at 1 AU is E = 250 g/cm 2 , which trans- 
lates to an initial disk mass of 2 x 10~ 3 Mi. The aspect 
ratio of the disk (h) is chosen to be 0.05. The temperature 



profile is T(r) cx r~ l , which is the result of the assumed 
constant aspect ratio and it is kept unchanged during the 
simulations (no energy equation is solved). 

Due to the mass loss during the initial stages of the sim- 
ulation, the disk mass is roughly 1.5 x 10~ 3 Mi, when the 
quasi-equilibrium is reached. Mass loss occurs in all sim- 
ulations where a secondary is present, however it is more 
severe for an eccentric secondary. If the secondary is on a 
circular orbit, the eccentricity of the gas disk is increased 
in the outer disk (see Fig. [TJd). The reason for mass loss 
is the following. Let us consider a fluid parcel that is ini- 
tially in a circular orbit in one of the outermost grids of the 
simulation. Due to the perturbation of the secondary, this 
fluid parcel is forced to an elliptic orbit. As the apoastron 
of the fluid parcel is now greater than its original orbital 
distance, the fluid parcel leaves the computational domain 
through the open boundary and mass loss occurs. The prob- 
lem is more pronounced for an eccentric secondary. Figure 
[2Jd shows the eccentricity of the gas disk when the sec- 
ondary is in apoastron. However, the gas eccentricity in 
(and shortly after) periastron is above unity, therefore the 
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Table 1. Initial parameters and results of the FARGO and dust simulations. 



ID 


Mi 


M 2 


o 


a. 


a 




r 


Inner b. 


t 


<£> 


< 77lin > 


< 777™,+ > 


< t s > 












[AU1 






fAUl 




[orbits] 


fsrl 

L&J 


L&J 






'single' 


1.4 











10" 


12 


0.4 


non-refl. 





0.004 


264 


446 


7.9 x 10" 


-■> 


'no_ecc' 


1.4 


0.4 





20 


10" 


2 


0.4 


non-refl. 


100 


0.012 


4.3 


1.8 x 10" 2 


2.0 x 10" 


-2 


'fiducial' 


1.4 


0.4 


0.4 


20 


10" 


2 


0.4 


non-refl. 


100 


0.087 


1.6 


4.8 x 10" 4 


4.6 x 10" 


-3 


'openb' 


1.4 


0.4 


0.4 


20 


10" 


2 


0.4 


open 


100 


0.045 


2.8 


4.0 x 10" 4 


5.5 x 10" 


-3 


'alpha' 


1.4 


0.4 


0.4 


20 


10" 


3 


0.4 


open 


200 


0.126 


4.5 


2.4 x 10" 3 


8.1 x 10" 


-3 


'rjn' 


1.4 


0.4 


0.4 


20 


10" 


2 


0.1 


non-refl. 


100 


0.087 


6.1 


6.4 x 10" 4 


5.2 x 10" 


-3 


'sep' 


1.4 


0.4 


0.4 


15 


10" 


2 


0.4 


non-refl. 


200 


0.059 


0.88 


7.7 x 10" 7 


2.4 x 10" 


-3 


'sec_m' 


1.4 


0.6 


0.4 


20 


10" 


2 


0.4 


non-refl. 


100 


0.093 


0.48 


3.8 x 10 -4 


1.8 x 10" 


-3 



In this table, Col. 1 describes the identifier of the simulation pointing out the parameter which is different compared to the 'fiducial' 
simulation (e.g., 'alpha' means that the turbulence parameter (a) is different than in the 'fiducial' simulation). Cols. 2 and 3 are 
the masses of the primary and secondary respectively, Cols. 4 and 5 are the eccentricity and semi major axis of the secondary 
respectively, Column 6 is the turbulence parameter (a), col. 7 is the radius of the inner boundary. Column 8 describes the boundary 
condition applied at the inner boundary, it can be non-reflective or open. Column 9 is the number of secondary orbits simulated 
before the equilibrium of the gas is reached, col. 10 is the surface density weighted average eccentricity of the gas at the end of 
the simulations. The properties of the final dust populations using the erosion model is described in cols. 11-13. Column 11 is the 
average mass of the particles between 0.5 and 1 AU, col. 12 is the average mass of particles between 4.5 and 5.5 AU, col. 13 is the 
average stopping time throughout the disk. 



gas in the outer disk is on a hyperbolic orbit and the gas 
disk between 6 and 8 AU is practically depleted (Fig. [5^). 

2.1.2. Boundary conditions 

We use most of the time the so called non-reflective bound- 
ary condition at rj„ which effici ently remove s any reflected 
waves from the inner boundary ([Crida et all 120081 ) . At ev- 
ery time step, the density in the zeroth ring is set to the den- 
sity of the first ring, which is rotated to simulate wave prop- 
agation and avoid wave reflection. The open inner boundary 
condition lets matter freely flow out from the grid. 

We always use an open boundary condition at r out , 
which lets the matter leave the computational domain dur- 
ing the initial stages of the simulation. This boundary con- 
dition sets the radial velocity to zero, if it points inwards. 
The radial velocity is extrapolated otherwise and the mat- 
ter can flow out. We also set a density floor in the outer 
ring, which we choose to be 10 orders of magnitude smaller 
than the average initial density. As explained in the pre- 
vious Section, the eccentricity in the outer disk is above 
1 when the secondary is in periastron. Therefore matter 
is continuously lost through the outer boundary and the 
density is decreasing after every passage of the secondary. 
We find that when the density is below a critical value 
(roughly 15-20 orders of magnitude below the average den- 
sity), negative density values occur due to numerical errors. 
The density floor criteria ensures that this does not happen 
and its effect on the surface density profile is undetectable. 

2.2. Quasi-equilibrium of the gas disk 

The structure of a disk around a single star is different 
from a disk which is subject to the perturbation of a close 
secondary, therefore we first follow the relaxation of the disk 
into a new quasi-equilibrium before modeling the motion of 
the particles in a binary system. 

Figures [TJi and b show the surface density profile and ec- 
centricity profile of the gas at different times in the 'no_ecc' 
simulation respectively. We see that the disk material re- 



structures, the outer radius of the disk becomes 7 AU after 
100 orbits. We do not see strong waves (spiral arms) in 
the disk, we see however that the eccentricity of the disk 
increases towards the outer part of the disk. 

The eccentric secondary excites two spiral waves when- 
ever it passes periastron. The waves propagate all the way 
to the inner disk. During the initial orbits of the secondary, 
the spiral arms carry a large amount of matter outside the 
computational domain which is then lost from the system 
due to the increased gas eccentricity at periastron. Between 
two periastrons the disk becomes more circular, it relaxes. 
The disk radius is truncated to approximately 5-6 AU as 
shown in Fig. [2ja, for the 'fiducial' simulation and the disk 
material restructures to a new equilibrium after 60 orbits. 
The azimuthally averaged eccentricity of the gas is pre- 
sented in Fig.[2}D. The average value is 0.09. The eccentricity 
increases between r = 4 and 8 AU due to the stronger per- 
turbations from the secondary. Table [T] shows the required 
number of orbits to simulate in order to reach equilibrium, 
the gas and the surface density weighted average eccentric- 
ity at the end of the simulations. The eccentricity of the 
gas is low, if the companion is on a circular orbit. The open 
inner boundary reduces the eccentricity by 50%, the lower 
a value increases the eccentricity by 50 %. The smaller in- 
ner boundary radius has no effect on the eccentricity. It is, 
however, interesting that a closer companion reduces the 
eccentricity of the gas compared to the 'fiducial' simula- 
tion. 

We note that although we mostly use a non-reflecting 
boundary condition, we do not favor this boundary condi- 
tion over any other possibility. The exact amount of mass 
deposited from the disk onto the star is determined by the 
interaction of ionized gas at several stellar radii from the 
star and the magnetic fields. It is far from clear what bound- 
ary condition should be used in a hydrodynamical simula- 
tion where the (sim ulated) inner disk ext ends to 0.1 AU 
only. As an example iKlev fc Nelson! (|2008l ) uses four differ- 
ent types of boundary conditions and concluded that the 
different boundary conditions result in different disk prop- 
erties. This is what we also observe. The built in boundary 
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conditions of FARGO are rigid, open and non-reflecting. 
The rigid boundary condition introduces reflected waves, 
the fully open boundary condition results probably in too 
high accretion rates. Therefore we prefer the non-reflecting 
boundary condition, but our choice is not based on more 
solid arguments than the above. A more thorough approach 
would be to explore the effects of all boundary conditions. 
However, we concentrate on the dust physics in this paper, 
we do not feel that such an approach is necessary. 

2.3. Dust 

In this Section we review the equations of motion of a 
dust particle, the numerical method used to solve these 
equations and we present the tests we performed to ver- 
ify the validity of our scheme. We also describe the erosion 
model used in our model. The full collision model, including 
growth, is described in Sec. 13.3.11 



2.3.1. Equations of motion 

A dust particle feels the gravitational force of the primary 
and secondary stars, the aerodynamical drag and inertial 
forces. The particles considered in this work have negligible 
masses compared to the masses of the binary stars, thus we 
can use the equations of the restricted three-body problem. 
The equation of motion of a dust particle in a spherical 
coordinate system (r, ip) is 



dr 
~dt 



~dt 



L/r 2 



dv r 
~~dt 



d<S> 

dr 



~t~ * d ,r ~t~ 



d --- d - + F d + F 



(3) 



(4) 



(5) 



(6) 



where v r is the radial velocity of the particle, L is the spe- 
cific angular momentum, $ is the common gravitational 
potential of the two stars (<I> = —GMi/r — GAl2/r s , where 
r s is the distance from the secondary star), Fd. r and Fd }tp 
are the radial and azimuthal drag forces respectively (see 
description in the next paragraph), and Fi n ^ r and Fi ntV are 
the radial and azimuthal inertial forces, respectively (see 
description in a later paragraph). 
The drag force is calculated as 



F d = At;, 



(7) 



IWeidenschillind (jl977t) ). We refer the reader to the above 
mentioned papers for further details of the calculation. 

The inertial forces arise since the coordinate system in 
FARGO is centered in the primary star and not in the cen- 
ter of mass of the system. Therefore, the center of mass 
orbits around the primary star. The inertial forces are cal- 
culated in Cartesian coordinates for simplicity and later on 
they are transformed into spherical coordinates. The x and 
y components of the inertial forces in such a system are 

F in , x = i? c cos$ c +2i? c $ c sin$ c -$^ c cos$ c -$ c i? c sin$ c , 

( 8 ) 

F in .y = R c sin $ c -2ii c $ c cos <S> C -<$> 2 C R C sin $ c +<£ c i? c cos $ c , 

(9) 

where R c and $ c are the (time dependent) distance and 
angle of the center of mass relative to the primary star. In 
case the secondary star is on a circular orbit (that is R c (t) = 
const, and $ c (i) = wi), the inertial forces are reduced to 
F- 



-lo R c cos <1> C and F;, 



-LJ i? c sin <& c , which 



have the same form as the centrifugal force expression. In 
the case of an elliptic secondary, the full expressions in Eqs. 
[8] and [9] must be used. 

2.3.2. Numerical method 

We use a sta ndard first orde r leapf rog integrator similar to 
that used by iPaardekooperl (|2007t ) . The discretized forms 
of Eqs. [SHU upon solving the drag force implicitly are 



V r0 



r = r + v r oAt 
tp = ip + L At/rl 



Atf2 fc 



L Q + At 



L = 



,rCt k 



Fan 



1 



Atn,, 



(10) 
(11) 

(12) 
(13) 



where At is the time step, u gir and the radial and 

azimuthal gas velocities respectively. Solving the drag force 
implicitly enables us to simulate the motion of an arbitrary 
small particle. We obtain the gas density and velocity val- 
ues from the FARGO code. The gas frames are saved 100 
times during one orbit of the companion, but the time step 
of our method is much smaller than that (usually 10~ 3 
yr). Therefore we linearly interpolate between two neigh- 
boring time frames (ti and t i+ i) to obtain the gas proper- 
ties at an intermediate time of the simulation t S i m , where 
U < isim < U+i- To determine the exact gas density and ve- 
locity components at the position of a particle, we perform 
bilinear interpolation using the four neighboring values. 



v g is the 



where flk is the Kepler frequency, Av = v d 
relative velocity between the dust and the gas and t s is 
the stopping time of the particle. The stopping time (or 
friction time) is the time the particle needs to react to 
the changes in the motion of the su rrounding gas. We fo l- 
low the de s criptio n used already by Paardekooperl <|2007t> . 
iLvra et all (|2009t ) or IWoitke & Hellind (|2003l ) to calculate 
t s . In this description, the stopping time for particles of in- 
termediate sizes is calculated as an interpolation between 
the Epstein regime (for par ticles smaller th an the mean free 
path of molecules in the gas lEpsteinl (|1924D ). and the Stokes 
regime (for particles much bigger than the mean free path 



2.3.3. Initial and boundary conditions 

The particles are distributed randomly in the disk in a way 
that the initial surface density profile of the dust is the 
same as the gas surface density profile. The particles have 
Kepler velocities initially. The boundary conditions for the 
dust is open both at the inner and at the outer boundary. If 
a particle leaves the FARGO grid, the particle is removed 
from the simulation. This differs from the gas boundary 
conditions which leads to some discrepancy between the 
surface density profile of the gas and small (t s ~ 0) particles 
(see Sect.^XP- 
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Fig. 3. The surface density of the gas (contour figures on the right sides) and the dust (left sides) initially at t = 100 
orbits (a) and after t = 105.6 orbits of the companion (b). The spiral waves in the gas as well as in the dust are clearly 
visible on Fig. b. 



2.3.4. Tests 

We performed three tests to ensure that our numerical 
method is correct. 

Stopping time of particles equals zero. In this case the parti- 
cles are well coupled to the gas. Thus the gas to dust ratio 
has to be constant during the simulation. We placed 10 5 
particles with t s = randomly in the hydro grid after the 
gas reached quasi-equilibrium (after t = 100 orbits). The 
initial distribution of particles is the same as of the initial 
gas distribution (see Fig. [3k.). We follow the motion of these 
tracer particles for 334 yrs (5 orbits) after which we com- 
pare the gas and dust distribution again (see Fig. |3Jd). The 
tracer particles concentrate in the spiral arms similarly to 
the gas, but some differences can be observed. The inner 
part of the dust disk is depleted compared to the gas disk 
due to the different boundary conditions used for the gas 
and the dust. Although the inner radius of the gas disk is 
at 0.4 AU (unless otherwise stated), we only consider dust 
particles outside 0.5 AU. 

The different boundary conditions for the gas and the 
dust is justified because the dust behaves like a pressure- 
less fluid. Although distant gas parcels can communicate 
with each other via pressure effects, there is no interaction 
between dust particles at different parts of the disk (i.e. the 
motion of any dust particle is determined by its local gas 
density and velocity only). 

Other differences can be explained by the presence of 
shocks in the gas. Whenever shocks appear (secondary in 
periastron), jumps in the radial and azimuthal velocity is 
present. As we use bilinear interpolation to determine the 
gas velocity at the position of the dust particle, these jumps 
(shocks) are smeared out, therefore the trajectory of the 
tracer particles somewhat differs from the trajectory of the 
gas. Higher grid resolution helps this problem, but as the 
shocks cannot be resolved in hydrodynamical simulations, 
some difference between the gas and t s = particles is 
always present. As we are not interested in the evolution of 
the dust surface density or the gas to dust ratio, this effect 
is not a serious problem in our investigation. 

Intermediate stopping time We compare th e drift speed of a 
test particle with the analytical results of iTakeuchi fc Linl 



(2002). If the gas density profile is constant over the radius 
and the temperature profile follows a T(r) oc r _1 , the radial 
drift velocity of the particle compared to the gas (Av r ) can 
be given as 



where h is the aspect ratio of the disk, h — 0.05 in our sim- 
ulation. Fig.@]shows the drift velocity values for t s = 0.005 
particles with the '+' signs, the analytical drift velocity dic- 
tated by Eq. [14] is the solid line. The difference between the 
analytical and calculated drift speeds is at least six orders 
of magnitude smaller than the drift speed itself, therefore 
we conclude that our implicit integration scheme is correct 
and the obtained relative velocities are sufficiently accurate. 

Infinite stopping time We placed a Jupiter mass planet ini- 
tially at (1,0) coordinates on a circular orbit (no secondary 
star). In case of a particle with t s 3> 1, the drag forces 
are negligible and the three body problem determines the 
motion of the particle. Thus a particle placed at the L4 or 
L5 Lagrange points should remain in the vicinity of these 
points and librate. Our test particle placed in the L4 point 
reproduces this movement. 

2.4. Erosion model 

As described in the previous sections, we follow the motion 
of individual large dust aggregates, measure their velocity 
relative to the gas. Naturally, we can only integrate the mo- 
tion of a limited number of particles. We assume that the 
rest of the dust population exist as "background" particles 
having small sizes (~ Ifim - small dust population). This 
small dust population is very well coupled to the gas (their 
relative velocity is practically zero) therefore it is reason- 
able to assume that the relative velocity we measure is the 
collision speed between the large simulated particles and 
the small dust population. Our collision model is based on 
this assumption. 

One could also obtain the relative velocity between 
two larger aggregates, but one has to face some problems. 
The distance between two aggregates is not negligible. The 
Kepler shear, which is present on even sub-grid scales, in- 
troduces relative velo cities which can e asily be several tens 
of meters per second (jLvra et all I2009F) . One has to remove 
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Fig. 4. Radial velocity difference between a dust particle 
with t s = 0.005 and the gas. The solid line is the analytical 
drift speed, the '+' signs are obtained using our integration 
scheme. 

all such effects to obtain true collision speeds (e.g., relative 
velocity at zero distance). We discuss these effects in more 
detail in Sect. 

Our collision model includes erosion, if the measured 
relative velocity is greater than 1 m/s. First, we calculate 
how many micron-sized monomer (from the small dust pop- 
ulation) collides with our aggregate during a time step. We 
calculate the aerodynamical cross section of the aggregate 
as 

a cr =irr 2 , (15) 

where r is the radius of the aggregate. We calculate the 
total gas mass that the aggregate swept through 



V x p gas = a cr AvAt x p gas , 



where V is the gas volume the aggregate travelled through, 
At is the time step of the calculation, and p gas is the gas 
density. We assume a typical gas to dust ratio of 100 for 
the small dust population. Although we follow the motion 
of some individual large particles, their mass compared to 
the total dust mass is negligible, therefore we can assume 
that all dust mass is present in the form of micron sized 
monomers. The number of collisions during a time step is 
then given by 



"-coll 



mo 



where mo is the mass of a monomer. Based on the lab- 
oratory experiments with erosion of du st aggregates, the 
mass of the particl e after a collision is (jGiittler et al.l . l2010t 
IWurm et al 



particle 



m i d 



8 Aw 

— Win 

60 lOOcm/s 



(18) 



Using this equation, however, resulted in very slow mass 
loss, thus the integration of hundreds of secondary orbits 
would have been necessary to reach the hnal particle sizes. 
We therefore speed up erosion artificially, as we are only 
interested in the final particle sizes and not in the exact 
erosion timescales. Our erosion equation is 



TO id - 10m 



Av 



We obtain final particle masses in less than 20 orbits using 
the equation above. 

This is clearly a simplified collision model as in a realis- 
tic dust distribution, the collision partners can have higher 
masses, thus the number of collisions per time step is re- 
duced. Also, the collision physics between similar sized par- 
ticles is different (more about this in Sec. 14. 2p . 

We note that after an erosion event, the cross section 
and the radius of the particle, and therefore the stopping 
time, the relative velocity and the number of collisions 
changes. If multiple erosion events happen during a sin- 
gle time step, these changes can be unresolved. To avoid 
such a numerical artifact, we use an adaptive time step- 
ping scheme which reduces the time step by a factor of 
two, when the mass-decrease of the particle (m new /Am) is 
greater than 10~ 3 . Also, the time step is increased by a fac- 
tor of two up to a maximum value of 10~ 3 , if m new /Am is 
greater than 10~ 3 . We find, however, that the time-step is 
never decreased during the simulations even when Eq. [T5] 
is used. 

In principle, a simulation would be completed when the 
mass of the dust particles do not change anymore. This 
condition could be fulfilled if the particles had a fixed or- 
bit. However, due to radial drift, this is not the case. By 
performing numerical experiments, we find that after 20 or- 
bits of the secondary, a mass distribution that only slowly 
varies due to radial drift is achieved, if the starting size 
(or mass) of the particles is not too far from the final dust 
mass. The initial dust mass therefore is not the same for 
all the simulations. In the 'single', 'no_ecc' and 'fiducial' 
simulations, the initial particle mass is 500 g, 100 g and 
50 g respectively. As seen on Figs. [SJ H2 and [JJ all particles 
lose mass by the end of the simulation, therefore the initial 
mass values are appropriate. 



(!6) 3. Results 



100 cm/s' 



(19) 



3.1. Dust in disks with and without a secondary 

We proceed in three steps in order to understand the effects 
of an eccentric binary on the dust population. First, we de- 
termine the properties of the dust particles in a disk around 
a single star ('single' set-up). Then, we compare these re- 
sults with a disk that is restructured by a non-eccentric 
secondary ('no_ecc' set-up), and finally, we examine the ef- 
fects of an eccentric binary ('fiducial' set-up). 

The gas density and velocity profile in the 'single' model 
is the initial disk profile used for all of the simulations. 
That is a power-law disk profile, which is in hydrostatic 
equilibrium and has a constant accretion-rate throughout 
the disk (for details see Sec. 12.1. ip . We keep this density 
and velocity profile constant during the simulation of the 
particles. As shown in Tab. [1] the eccentricity of the gas 
is the lowest of all set-ups (0.004). The small eccentricity 
is the result of the radial drift motion of the gas and it is 
also constant in time. The only source of relative velocity 
in this model is the differential radial drift of the dust. 
The dust simulation is started with particles having 500 
g. We measure time in units of the orbital period of the 
secondary in the 'no_ecc' or 'fiducial' models to be able 
to easily compare the results of the different models. The 
mass of the particles is kept constant for one orbital period 
(66.7 yrs) during which time the particles forget their initial 
Kepler velocity and sufficiently couple to the gas. After 66.7 
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Fig. 5. The evolution of the mass of the particles in case 
of the 'single' set up (no secondary present, the disk is not 
perturbed). The x axis is time in orbital periods of the 
secondary in the 'no_ecc' or 'fiducial' set ups (one orbit 
takes 66.7 yrs). The y axis is the mass of the particles. The 
colors represent how many particles are present at a given 
location of the parameter space. 
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Fig. 6. Same as Fig. [5] for the 'no_ecc' set up (a secondary 
is present on a circular orbit). The disk is truncated, which 
results in lower particle masses compared to the 'single' 
simulation. 



yrs we turn on erosion and follow the mass-change of the 
particles as shown in Fig \b\ . Figures [8^l and b show the 
final mass and stopping time of the particles as a function 
of their distance from the primary star respectively. The 
stopping time in general does not change strongly with the 
distance, therefore it can be characterized with the average. 
However, the mass of the particles can vary greatly with 
distance, thus we characterize the mass with its average 
around the inner and the outer edge of the disk. The average 
mass between 0.5 and 1 AU as well as between 4.5 and 5.5 
AU, and average stopping time of the particles are shown 
in Tab. [Q 

A similar procedure is followed for the 'no_ecc' simula- 
tion, where a secondary is on a circular orbit. The disk is 
truncated, the outer edge of the disk is around 7 AU. The 
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Fig. 7. Same as Fig. [5] for the 'fiducial' simulation, where 
a secondary is present on an eccentric orbit with e = 0.4. 
The secondary closely approaches the disk at periastron (at 
every half orbital time) at which times the particles at the 
outer disk experience increased relative velocities, therefore 
they erode. 



relative velocity in this, and all the other simulations is the 
result of both the radial drift and the ever-changing gas 
velocity profile. The particles have initially 10 g mass. The 
mass-evolution of the particles in this simulation is shown 
in Fig. [6] The final masses and stopping times of the par- 
ticles as the function of their distance are shown in Figs. 
[5J; and d respectively. The average mass at the inner, and 
outer disk as well as the average stopping time can be seen 
in Tab. HJ 

The mean free path of the gas atoms and molecules is 
of the order of 1 meter in these disk models due to the 
moderate densities. Furthermore the radius of the particles 
is below 10 centimeters. Therefore, we conclude that the 
particles are in the Epstein regime. The stopping time in 
this regime is 

3w ,„ . 

t s =t Ep = - , (20 

■IVthPgVcr 

where v t h is the thermal velocity of the gas. The relative 
velocity of the particle depends on how fast the particle 
can adapt to the changes in the gas velocity field, e.g., the 
relative velocity is higher, if the stopping time is higher. 
However, the gas density is decreased in a truncated disk. 
As seen in Eq. [301 lower gas density increases the stopping 
time, thus the relative velocity is higher as well. Erosion, 
therefore, reduces the masses of the particles in the outer 
disk. 

Therefore, one effect of the secondary is the reduced 
particle masses in the outer disk. The stopping time of the 
particles is also lower than in the 'single' model. As seen in 
Tab. HI the eccentricity of the gas is three times higher in 
this simulation, therefore the particles, which would move 
on a circular orbit, experience stronger changes in the gas 
velocity field. Thus the stopping time is also reduced until 
the relative velocity reaches the critical 1 m/s value. The 
destructing effect of the gas eccentricity on the dust popu- 
lation affects the entire disk. 

The particles in the 'fiducial' simulation have initially 5 
g masses. The mass evolution is the most dramatic in this 
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case. The companion approaches the primary as close as 12 
AU in periastron, therefore the outer edge of the disk ex- 
periences strong disturbances periodically. Every time this 
happens, the mass of the particles in the outer disk is re- 
duced, which can be seen in Fig. [7J The secondary in pe- 
riastron excites spiral arms in the gas disk as it is pushed 
inward. Therefore, a companion in an eccentric orbit fur- 
ther increases the eccentricity of the gas, thus decreasing 
the mass and the stopping time of the particles (see Tab.[U 
Figs. [5^ and f). This effect influences the dust distribution 
in the whole disk. 

We conclude that the average stopping time is reduced 
by one order of magnitude in the 'fiducial' simulation com- 
pared to the 'single' simulation. The average mass of the 
particles however is roughly four orders of magnitude lower. 

3.2. Exploring the parameter space 

We perform altogether eight simulations with different ini- 
tial conditions as shown in Tab.[TJ We see that if the turbu- 
lent parameter is kept constant, different parameters (like 
the inner boundary condition, separation of the binary sys- 
tem, mass of the secondary) affect the gas, especially the 
eccentricity of the gas as shown in Tab. [T] We perform sim- 
ilar dust simulations as the ones discussed in the previous 
section in all of these models and measure the average dust 
particle mass at the inner and outer disk, as well as the 
average stopping time of the particles at the end of the 
simulation. The results are shown in Tab. [T] As a general 
trend, we can say that given a turbulence parameter, the 
average mass and stopping time of the particles decrease as 
the eccentricity of the gas increases. 

The 'alpha' simulation with ten times lower a value pro- 
duces the most eccentric disk, but the mass and stopping 
time of the particles are higher than the average of all the 
a = 10 -2 cases. The reason for this is that if the turbulence 
is lower, the radial flow of the gas is lower, the fluid behaves 
in a less viscous way. Therefore the dust relative velocities 
are also decreased allowing heavier particles to exist in such 
a disk. 

The radius of the inner disk in the 'r_in' simulation 
is close to the usually assumed sublimation limit of the 
dust (where the temperature is higher than 1500 K). We 
performed this special simulation to see whether the ef- 
fects of the secondary reaches the very inner parts of the 
disk. Previous work done on the growth of planetesimals 
in binary systems showed that the planetesimals in the in- 
ner regions are less affected by the secondary star, there- 
fore these regions are favorable for their growth ( see e.g., 
iThebault et al.l (|2004l )). In contrast, the growth of dust ag- 
gregates is less responsive to the direct gravitational per- 
turbation of the secondary, but reacts sensitively to the gas 
dynamics. As the eccentricity of the gas is pumped up by 
the secondary throughout the disk, even the close regions 
to the primary star do not favor dust growth. 

3.3. Long term prediction 

The erosion model alone can tell us the maximum particle 
sizes in the disk that does not experience any fragmenta- 
tion anymore. In the case of the 'fiducial' set-up, however, 
the relative velocity of the particles is only increased when 
the secondary is in periastron as shown in Fig. |H1 The rela- 
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Fig. 9. The relative velocity distribution of aggregates as a 
function of time in the 'fiducial' simulation. The x axis is 
in orbital time units, the y axis is the relative velocity in 
cm/s. The colors indicate how many particles are present 
at a given location of the parameter space. The white solid 
line shows the average relative velocity of all particles as a 
function of time. 

tive velocity between two periastrons is much smaller than 
the critical velocity. The particle could experience growth 
during these periods, which alters the equilibrium mass dis- 
tribution of the particles. 

To examine the long term evolution of the dust and 
calculate the equilibrium mass distribution as the function 
of the distance from the primary, we add coagulation to 
the collision model, distribute the particle masses randomly 
and follow the mass evolution of the particles for one orbit 
of the secondary. For more than one orbit, the radial drift of 
the large particles can be significant, thus non-local effects 
(due to radial drift) could alter the equilibrium mass distri- 
bution. Simulating less than one orbit is not sufficient since 
the appearance and dissipation of the spiral arms happens 
in an orbital period. We then calculate the growth/erosion 
rate of the particles {{m,2 — mi)/mi/T, where m\ is the 
starting mass, m>2 is the final mass of the particle in one 
orbit - T) as a function of their distance and stopping 
time and identify the zero erosion/ growth rate curve. If 
we assume that the gas density and velocity profiles do not 
change significantly over multiple orbits, this curve identi- 
fies the equilibrium mass distribution between growth and 
erosion. In the following section, we describe the changes 
in the collision model, and the new initial conditions. 

3.3.1. The complete collision model 

In this section we introduce a new collision model that in- 
corporates growth and erosion. Although we call this col- 
lision model complete, it does not incorporate all known 
collision types (see Sec. 14.21) . 

Growth Equation [T7J describes the number of collisions 
which takes place during one time step of the calcula- 
tio n. According t o the experiment-based collision model 
of iGuttler et ail IpOlCh . it is realistic to assume that a 
monomer sticks to an aggregate of arbitrary mass, if the 
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Fig. 8. The figures on the left side (a, c, e) show the mass of the particles as a function of their distance from the primary 
at t = 120 orbits for the 'single', 'no_ecc', and 'fiducial' simulations respectively. The figures on the right side (b, d, f) 
show the corresponding stopping times of the particles at the same time. 



relative velocity is lower than the critical fragmentation ve- order to obtain the final masses of the aggregates in a rea- 

locity. Therefore, in the new collision model, n co ii number sonable computational time. This is not necessary anymore, 

of monomers are added to the aggregate, if Aw is less than as we follow the mass evolution of particles only for one 
1 m/s. 



Erosion In the previous collision model (see Sec. 
artificially increased the erosion efficiency (see Eq. [1 



we 
in 



secondary orbit (see Sect. I3.3| ). Therefore we use Eg. [TSlto 
follow the erosion of the particles, if the collision velocity is 
greater than 1 m/s. 
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3.3.2. New initial conditions 

In the following calculations the particles are initially uni- 
formly distributed in radius between Ti n and r out (previ- 
ously, the particles were distributed according to the gas 
surface density profile). Also, their masses are initially dis- 
tributed uniformly in log-space between a monomer mass 
and 500 g. The initial velocity of the grains is Keplerian, 
and (as the stopping time of all particles is lower than the 
orbital time) for one obit we follow only the motion of par- 
ticles to let the particles couple sufficiently to the gas. We 
turn on the collision model after one orbit of the secondary. 

3.3.3. Results 

Figures ITOb and b, ITOb and d, and fTUb and f show the 
growth/erosion rates for the 'single', 'no_ecc', and 'fidu- 
cial' simulations, respectively. The growth/erosion rates 
are expressed in normalized dimensionless units ((to2 — 
ml)/mi/T, where mi is the initial mass of the particle, 
rri2 is the final mass, and T is the orbital period of the sec- 
ondary). If e.g., the growth rate is unity, it means that the 
large aggregate doubled its mass during one orbital time of 
the secondary. 

Figure [TUb and b for the 'single' set-up shows that the 
growth/erosion rate decreases/increases with radius which 
is due to the decreasing density, thus lower number of col- 
lisions, as a function of distance. 

Figures [TUb and d, and [TUb and f show more complex 
features. The growth/erosion rates first decrease/increase 
as a function of distance, but it increases/decreases again at 
the outer parts of the disks. This behavior can be explained 
by the eccentricity of the gas. As seen in Figs. [JJj and [2b, 
the eccentricity of the gas is higher in the outer parts of the 
disks due to the perturbation of the secondary. Due to the 
the eccentric gas at the outer disk, the relative velocities 
rise, therefore as long as the relative velocity is not above 
1 m/s, the growth rate is increased again. 

We compare the equilibrium curve on these figures with 
the obtained solutions in Figs. [5b-f. The agreement is good 
for the 'single' simulations. However, differences can be ob- 
served in the case of the 'no_ecc', and 'fiducial' simulation 
and in all other simulation set-ups where the secondary 
is eccentric. The stopping time and mass in the updated 
collision model is generally higher in the outer disk than 
in the previous collision model including only erosion. The 
spiral arms excited by the secondary erodes the particles 
in periastron, but the aggregates can grow significantly be- 
tween two periastrons. Therefore, calculating only erosion 
can underestimate the equilibrium size/stopping time of the 
particles. 

Other difference can be observed in the inner disks. 
The mass and stopping time of the particles are higher in 
the full collision model than in the previous (erosion only) 
model. There are two reasons for that. First, the particles 
in the previous model do not grow between two periastrons. 
Secondly, the radial drift of the particles become important 
so close to the inner disk edge. As we only follow the particle 
motion for one orbit in the full collision model, the drift of 
the particles is negligible. However, in the previous model, 
the motion and dust evolution of particles is followed for 20 
orbits (roughly 1500 yrs), therefore many particles in the 
inner disk with masses higher than 10 g were lost during 



the simulation. We give a more detailed discussion about 
the relevant timescales of the results in Sec. 14.11 

4. Discussion 

4.1. Validity of the results: relevant timescales 

Coagulation timescale We use such a collision model in Sec. 
13.11 where we start the simulation with centimeter-sized 
particles and we follow how they grind down to smaller 
particles. The question arises, in what timescale can these 
particles form by coagulation? Is it valid to assume that 
these particles are present in the disk already after 100 or- 
bits of the secondary? In order to answer this question, we 
estimate the coagulation timescale: 

^coag — 7 7 (21) 
n d Ava cr 

where rid is the number density of the dust aggregates with 
a given size (n d = p g /100/m, where p g is the gas density, 
m is the particle mass and we assumed a 100:1 gas to dust 
ratio.). Using typical density values at 1 AU in our models, 
assuming a particle mass of 1 g, and relative velocity of 10 
cm/s, we get that t coag = 3 x 10 3 yrs (50 orbits). Thus we 
conclude that particles can reach the initial particle sizes 
that we use in Sec l3.ll in 100 orbits, of course under ideal 
conditions. 

Viscous timescale of the gas The viscous timescale of the 
disk at a given distance is given by 

t vis = r 2 /v, (22) 

where r is the distance from the central star, v is the vis- 
cosity which is defined in Eq. [TJ The viscous timescale in 
our model at 1 AU distance is of the order of 10 4 yrs us- 
ing a = 10~ 2 , and 10 5 yrs assuming a = 10~ 3 . This cor- 
responds to 150 to 1500 secondary orbits depending on a. 
The viscous spreading and accretion of the disk restructures 
the density profile of the disk during one viscous timescale 
and our results have to be treated with caution if used for 
a longer timescale than this. 

Radial drift timescale of the dust Radial drift (vd) has two 
sources: the drift of individual particles with respect to the 
gas (vd) and drift caused by accretion processes of the gas 
{vd g ), thus the total radial drift velocity is vd — Vd + Vd g - 
We measure the total radial drift speed of particles with 
different Stokes numbers (or dimensionless stopping times) 
at 1 AU in the 'single' simulation. The radial drift speed 
measured in that model is not affected by the eccentricity 
variations of the gas disk. The radial drift timescale is then 
the time the particle needs to drift into the star, e.g., drift 
1 AU distance. It turns out that idrift = 2 x 10 4 yrs (300 
orbits), if the stopping time t s — 10~ 3 T p , where T p is the 
orbital period of the particle; idrift = 10 4 yrs (150 orbits), 
if t s = 10- 2 T p ; and i dr ift = 2 x 10 3 yrs (only 30 orbits!), if 
t s = 10 _1 T p . The particles located at the inner part of the 
disk in our simulations are lost into the primary within 150 
orbits. Although the radial drift is also a serious problem 
for single disks, this poses an even stronger constraint on 
planet formation in binary disks. As the disk in this binary 
system is small (the size limited by the separation of the 
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Fig. 10. The figures on the left side (a, c, e) show the growth/erosion rates of particles in the 'single', 'no_ecc' and 
'fiducial' simulations respectively. The x axis is the distance from the primary in AU units, the y axis is the mass of the 
particles in gram. The red contours (below the white line) illustrate the growth rate, white meaning that the particle 
doubled in mass during one orbit. The blue contours (above the white line) represents erosion rates, white in this case 
meaning that the particle completely eroded in one orbit. The figures on the right side (b, d, f) show the same quantities, 
but the stopping time is plotted on the y axis. The white solid line indicates the equilibrium between growth and erosion. 



system), there is no reservoir of dust mass available as in 
disks around single stars. Unless there is a mechanism to 
recycle the dust, or stop the drift of the dust before it gets 
accreted onto the star, the total dust mass of the disk is 
lost within several hundred orbits. Therefore, the available 
time to form planetesimals from the dust aggregates can be 
limited to a very narrow window. 



4.2. A critical look at the collision model 

The considered sources of relative velocity between the dust 
and the gas are the radial drift and the perturbations in the 
gas velocity field due to the secondary. In reality, however, 
particles are also stirred up by turbulent eddies, particles 
settle toward the midplane of the disk and these processes 
contribute significantly to the relative velocity between par- 
ticles. Therefore, our models underestimate the true rel- 
ative velocities by neglecting these effects, thus the final 
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particle masses and stopping times are probably somewhat 
overestimated. 

We do not consider collisions between similar sized par- 
ticles. Determining a realistic collision speed (e.g., relative 
velocity at distance = 0) of similar sized aggregates in hy- 
drodynamical simulations is non-trivial. One has to sub- 
tract the Kepler velocity from the velocity vectors of the 
two parti cles to get r id of the artifact of the Kepler shear 
(see e.g., iLvra et all ([2009)). Even after this, the relative 
velocity of the particles inside the same grid cell strongly 
depends on the distance between these aggregates. A true 
collision speed can only be determined by extrapolating the 
relative velocity to zero distances using sufficiently large 
number of particle pairs in a statistical way. We therefore 
only consider the relative velocity between tiny micron sized 
particles (that essentially move with the gas) and the dust. 
These relative velocities are defined at the exact same lo- 
cation, thus the velocities used throughout the paper are 
true collision speeds. 

Furthermore, we only consider two collision types: 
sweeping up monomers, and erosion by monomers. There 
are however nin e different coll i sion ty pes defined in the col- 
lision model of iGiittler et al.l (|2010[ ). The most critical of 
them is bouncing. As we only calculate collision speeds 
between a monomer and an aggregate, currently we are 
unable to consider any other collision type but the ones 
described above. However, collisions between similar sized 
bodies can be more destructive. Particles could fragment 
into smaller pieces, or the growth could halt due to bounc- 
ing ()Zsom et all . 12010ft . 

Considering all these limitations, our results can only be 
regarded as upper estimates of the final mass and stopping 
time distribution of dust in binary systems. 

4.3. Possible planetesimal formation processes 

We investigated dust coagulation/fragmentation in this 
work and found that the maximum aggregate sized formed 
by coagulation is roughly 1 g in the inner disk and 10~ 4 g 
in the outer disk in binary systems. We discuss in this sec- 
tion the possible methods to form planetesimals from these 
aggregates. We consider favorable places for coagulation in 
the disk, and also particle concentration mechanisms, where 
planetesimals can directly form from dust aggregates. 

Coagulation in a pressure bump. iKretke &: Linl (|2007l) ar- 
gued that a pressure bump forms around the snow-line, 
which stops the radial drift of the particles and concen- 
trates them. Planetesimal formatio n can proceed in suc h 
locations as shown by the work of iBrauer et al.l (|2008bft . 
If such a pressure bump exists in disks around primary 
stars, this is a favorable mech anism to form planetesi - 
mals. However, the recent work of lDzvurkevich et al.l (|2010t ) 
showed t hat the pressu r e bum p may be milder than pre- 
dicted by IKretke fc Lml ()2007t) . 

We performed two simulation with FARGO in which we 
simulated the stability of a simplified pressure bump in the 
presence of a circular and eccentric secondary. The initial 
surface density profile is altered to include a gauss-shaped 
'gap' (see[TTJ black line at t = orbits). The gas velocity 
is altered in a way that the surface density profile is stable 
on a timescale shorter than the viscous timescale. In order 
to make the pressure bump stable on any timescale, the 



kinematic viscosity of the gas has to be radial dependent, 
which is currently not possible with FARGO. We choose the 
a parameter to be a = 10~ 5 , thus the viscous timescale is 
~ 10 5 secondary orbits. The azimuthal velocity on the right 
side of the gap is greater than the Kepler velocity, therefore 
this pressure bump would be strong enough to stop the 
radial drift of the particles, thus capture them there. We 
performed a test simulation where no secondary was present 
and the pressure bump did not change significantly over 150 
orbits. 

Figure fTTh shows the evolution of the surface density in 
the presence of a secondary on a circular orbit. The sec- 
ondary pushes the matter inward and excites waves in the 
gas. Therefore, the gap is gradually filled with gas on a 
timescale which is much shorter than the viscous timescale 
of the disk. This picture is even more dramatic in case the 
secondary is on an eccentric orbit with e = 0.4 (see Fig. 
II li b) . The pressure bump disappears completely in 30 or- 
bits. Despite of the simplicity of these simulations, it shows 
that the pressure bump can evolve on a shorter timescale 
than the viscous or coagulation timescale, namely on the 
orbital timescale of the secondary. If the pressure bump is 
not stable for even a coagulation timescale, planetesimals 
cannot be formed in such a pressure bump. 

Other particle concentration mechanisms. The other 
planetesimal f o rmati on mechanism was outlined by 
iJohansen et a l. (2007). This scenario assumes that a large 
amount of the solid material is presented in dm-sized 
boulders (t s > 0.1) at the midplane of the disk. These 
boulders then concentrate in long-lived high pressure 
regions in the turbulent gas and initial over-densities 
are amplified further by the streaming instability. This 
mechanism forms 100 km sized objects on a very short 
timescale, however it is still debated whether this scenario 
wo rks for small e r part icles with t s < 10~ 2 . 

ICuzzi et al.l (|2008ft outlined an alternative concentra- 
tion mechanism to obtain gravitationally unstable clumps 
of particles, which can then undergo sedimentation and 
form a 'sandpile' planetesimal. In this model, turbulence 
causes dense concentrations of aerodynamically size-sorted, 
chondrule-sized particles. This mechanism effectively con- 
centrates even t s = 10~ 4 particles, therefore it seems a good 
candidate for planetesimal formation. However, this mech- 
anism requires a critical amount of dust over-density in the 
clumps. The periodic perturbation from the secondary stirs 
up the dust particles, therefore it is necessary to examine 
whether particles can concentrate sufficiently in spite of the 
sti rring of the secondar y. 

iBarge fc Sommerial (119951) : iKlahr fc Hennind (|1997[) 
and recently ILvra et al.l (|2008h showed that vortices accu- 
mulate dust particles very efficiently. Such vortices can be 
formed for example at the edge of the dead zone in disks. 
As discussed for the other concentration mechanism, it is 
necessary to examine whether the perturbation of the gas 
and the dust caused by the secondary does not prevent the 
concentration mechanism. 

5. Conclusions 

We have studied for the first time the coagula- 
tion/fragmentation processes of dust aggregates in a disk 
in binary systems. We have measured the relative velocity 
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Fig. 11. The surface density profile of a disk with an artificial pressure bump at t=0, 10, 30, 60, and 100 orbits in case 
of a secondary on a circular orbit (a) and on an eccentric orbit with e=0.4 (b). 



between the dust aggregates and the gas and assumed that 
this velocity is the same as the relative velocity between 
the dust and well-coupled monomers. Based on this, we 
have constructed a simple collision model. If the relative 
velocity is smaller than 1 m/s, the critical fragmentation 
velocity, the particle sweeps up monomers from the disk. If 
the relative velocity is greater than 1 m/s, the aggregate 
is eroded based on a recipe measured in laboratory exper- 
iments. Using this collision model, we have found that the 
eccentric secondary has two major impacts on the dust pop- 
ulation in the disk. 

— A secondary on a circular orbit truncates the disk, there- 
fore the particle masses/stopping times are decreased at 
the outer part of the disk. 

— A secondary increases the eccentricity of the gas disk. 
If the eccentricity of the gas (therefore the deviation 
of the gas velocity field from the equilibrium velocity 
field in accretion disks) is increased, the relative velocity 
between the dust and the gas is greater. This in turn 
results in even smaller aggregate sizes throughout the 
disk. 

Due to these effects, the average stopping time/mass of the 
particles in a disk with an eccentric secondary is reduced 
by one/four order/s of magnitude compared to the dust 
population which resides in a disk around a single star. 

We have also studied the possibilities to form planetes- 
imals from the dust aggregates investigates in this work 
(enhanced coagulation in a pressure bump and particle con- 
centration mechanisms). We found that a pressure bump 
can evolve on shorter timescales than the viscous or coag- 
ulation timescale, namely on the orbital timescale of the 
companion. This effect, if it is confirmed by more detailed 
hydrodynamical studies, may suggest that a pressure bump 
in a binary system is not long lived enough for the forma- 
tion of planetesimals. Particle concentration mechanisms 
in general need to cope with the continuous gravitational 
stirring and perturbations of the secondary. 
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